Jto 


_  >, 
E  a 
o  o 

■Q  * 
8  W 

■3? 

2  < 
a 

a>  w 
X.  o 
CO 


^  j)  £  5^  /  3-^ 


TITLE*  Theory  of  Human  Vibration  Response 


FRED  PRADKO,  RICHARD  LEE— aBjUilCJ 
U.S.  Army  Tank-Automoti  v| 

Warren,  Michigan  48090  —  — —  * 


?  vv»g?  -  «  | 

8  8  VF*tl 


s"r,f“n 

t—  to*  •*■  fc.'  m 


ABSTRACT  I 


LK& 


LIBRARY 
:  COPY 


Analytical  and  experimental  studies  of  whole  body  human 
dynamics  under  random  vibration  are  presented.  A  previously 
un reported  criteria  "absorbed  power11  is  developed  through  the 
application  of  transfer  functions.  This  method  is  applicable 
to  assess  the  effects  of  stationary  and  non -station ary  vibration 
records.  It  is  equally  effective  for  synthesis  or  analysis. 

The  linearity  of  human  response  to  vibration  is  established  on 
a  qualitative  basis.  . 

Examination  of  the  "absorbed  power"  criteria  indicates 
advantages  not  present  in  the  acceleration  measurement. 

"Absorbed  power"  is  a  scala**  quantity  which  may  be  described  by 
magnitude  only.  It  is  additive  and  may  be  summed  in  multi - 
degree-of-freedom  environments. 

"Absorbed  power"  does  not  require  frequency  spectrum 
analysis.  Optimization  studies  involving  human  dynamics  may  be 
conducted  continuously  without  the  time  lapse  that  occurs  for 
frequency  spectrum  analysis.  The  findings  of  this  research 
indicate  with  substantial  credence  that  "absorbed  power"  and  transfer 
function  techniques  may  provide  the  basis  for  a  universally 
usable  human  vibration  measurement  method  which  shal 1  be  applicable 
to  air,  sea,  and  land  transportation  media. 
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THEORY  OF  HUMAN  VIBRATION  RESPONSE 


FRED  PRADKO,  RICHARD  LEE,  and  VICTOR  KALUZA 
U.S.  ARMY  TANK-AUTOMOTIVE  CENTER 
WARREN,  MICHIGAN 

Int roducti on :  The  subject  of  human  vibration  has 
been  treated  by  many  investigators  with  some  of  the  earliest  work 
dating  back  to  1818.  Current  bibliographies  indicate  over  1500 
separate  studies  have  been  directed  to  the  problem  of  identifying 
human  response  to  vibration. 

In  1962,  it  appeared  that  existing  research  would 
continue  to  place  heavy  emphasis  upon  experimental  tests  to  assess 
the  effects  of  mechanical  vibration.  New  programs  based  upon 
analytical  approaches  were  not  appearing.  Therefore,  in  January 
19^3,  due  to  an  extreme  internal  need  for  forecasting  human 
response  to  vehicle  vibration,  the  U.S.  Army  Tank-Automotive 
Center  began  a  basic  research' program. 
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The  premise  of  this  work  was  oriented  towards 
random  vibration,  as  it  has  been  generally  recognized  that  random 
vibration  best  describes  the  environment  created  by  military,  air, 
sea,  and  automotive  transportati  on. 

This  paper  discusses  whole  body  human  response  to 
mechanical  vibration  below  60  cps.  The  theoretical  considerations 
are  based  on  the  view  that  man's  response  in  a  vibratory  environ¬ 
ment  can  be  determined  through  measurement  of  input  conditions 
only.  Separation  of  the  total  problem  (anatomical,  psychological, 
and  physiological)  into  this  narrow  premise  has  produced  a  new 
parameter  identified  as  "absorbed  power". 

Tjiis  parameter  is  completely  different.  It  has 
been  developed  to  describe  human  d/nami  cs  without  prior  knowledge, 
of  the  frequency  spectrum  and  can  identify  human  response  with 
single  numerical  distinctness. 


Average  "absorbed  power"  =  Lim 
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F(t)V(t)dt 
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Pi scussiont  The  word  "random"  is  used  primarily 
in  a  statistical  sense  and  implies  vibration  which  does  not  show 
frequency  or  amplitude  order  except  as  a  statistic.  Random 
vibration  in  general  describes  envi ronmental  excitations  which 
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exist  over  a  band  of  frequencies.  When  motion  is  periodic,  the 
accelerations  are  determinable  by  historical  examination.  When 
the  motion  is  random,  non-periodic,  only  the  probability  of  future 
-^occurrence  can  be  stated. 

Previous  to  this  research,  to  conveniently  cope  with 
random  vibration  calculations  and  experimental  work,  the  concept 
known  as  acceleration  density  was  used.  In  brief  review,  this 
concept  takes  the  square  of  the  RMS  value  of  acceleration  over 
specific  frequency  intervals  and  yields  the  mean  square  value  of 
the  acceleration  for  that  particular  bandwi th.  If  this  mean  squared 
value  of  acceleration  is  divided  by  the  bandwidth  and  the  bandwi th 
is  made  to  approach  zero,  the  resulting  value  can  be  plotted  versus 
frequency  over  a  spectrum  of  interest.  This  is  known  as  the 
POWER  SPECTRAL  DENSITY  (PSO)V  The  mechanical  analogy  of  pqwer  is  the 
rate  of  doing  work.  PSD  is  proportional  to  the  square  of  the 
vibration  amplitude.  Hence,  a  plot  of  PSD  versus  frequency  shows 
the  power  distribution  of  the  vibration  environment.  PSD  and  "ABSORBED 
POWER"  will  again  be  frequently  referred  to.  These  two  quantities 
are  not  interchangeable.  Meaning  and  application  of  these  two  terms 
are  distinctly  different.  PSD  statistically  describes  the  vibration. 
"Absorbed  Power"  measures  the  rate  at  which  the  vibratory  input  is 
used.  Normally,  comprehensive  descriptions  of  random  vibration  tests 
requi  re  PSD  specifications.  For  "white  noise"  random  vibration,  the 
PSD  is  flat;  that  is,  all  frequency  components  are  present  at  equal 
energy  levels.^ 


Power  spectral  density  and  other  commonly  used 
techniques  such  as  time  histories,  peak  distribution,  or 
probability  records  describe  the  vibratory  input  to  the  "occupant"  o 
of  the  "man"  but  fail  to  identify  his  subjective  response  to  this 
envi ronment. 

Man  has  the  capacity  of  evaluating  situations 
involving  touch,  taste,  sight,  and  sound.  Other  research  has 
successfully  established  the  property  of  "quantity"  in  many  of  these 
common  sensations,  but  the  problem  of  establishing  a  unit  to 
quantify  and  measure  the  perception  of  vibration  analytically  has 
not  been  previously  disclosed. 

Approach?  A  principal  concept  in  modem  analysis 
of  shock  and  vibration  is  the  transfer  function.  This  function  is 
the  complex  ratio  of  output  information  to  the  input  which  caused 
the  output. 

The  transfer  function  is  usually  denoted  G(S)  where 
S  =  q*  +  ju  in  the  Laplace  transform  of  a  function  g(t).  Sigma  (g-) 
represents  the  transient  term  and  jto,  the  steady  state  condition. 

The  transform  is  denoted  by  the  symbol  L[g(t)]  or  G($)  and  is 
defined  by  the  integral. 


•  co 


L[g(t)]  =  G(S)  = 


eStg(t)dt, 


v 

General  form  of  the  transfer  function  G(S)  is  the  ratio, 

nt  e \  '  _  Response  Function  ,  . 

Input  Function  '  ' 
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Each  combination  of  output  and  input  in  a  linear 
system  may  be  described  by  a  transfer  function.  The  transfer 
function  relating  the  acceleration  input  A.  to  the  output  A  is 
-^given  by  Equation  4.  The  transfer  function  relating  the  force  of 
vibration  F . n( S )  to  the  output  acceleration  is  given  by  Equation  5. 


Us)  «V“£<s>.. 

in 


(4) 


Gets)  =  (S). 

out 


(5) 


The  first  hypothesis  of  this  program  considered  the 
application  of  transfer  function  techniques  to  human  dynamics. 

The  essential  constraint  controlling  the  use  of  the 
transfer  function  is  that  the  system  in  question  responds  linearily 
within  the  bounds  of  interest.  Therefore,  application  of  the 
transfer  function  technique  to  human  dynamics  required  the 
establishment  or  acceptance  that  man  behaves  as  a  linear  system. 

The  frequency  spectrum  is  used  to  produce  the  transfer 
function.  In  this  research  program,  a  common  procedure  was  used; 


predetermined  sinusoidal  inputs  of  specified  magnitude  and  frequency 
were  produced  and  the  magnitude  and  phase  of  the  steady-state  output 
recorded. 

The  necessary  data  for  the  transfer  function  was 
collected  in  a  series  of  experiments  for  sinusoidal  and  random 
environments  involving  thirty-one  male  test  subjects  representing 
enlisted  military  personnel  from  the  U.S.  8th  Artillery  Group, 
Detroit,  Michigan,  and  civilian  personnel  from  the'  U.S.  Army  Tank- 
Automotive  Center.  The  statistical  summary  of  the  test  population 
is  provided  in  Appendix  A. 

The  data  collected  described  the  input  to  the  man 
in  terms  of  force,  acceleration,  and  frequency.  The  output  was 
measured  and  recorded  at  the  subject's  head. 

Figure  1  describes  the  vertical  mode  data  as  the  two- 
sided  interval  which  can  be  expected  to  bracket  the  true  mean  with 
90%  confidence.  The  ordinate  describes  the  RMS  (root  mean  square) 
value  of  acceleration  in  g's  for  1  through  30  CPS  under  the  following 
test  controls.  The  frequency  was  pre-set  and  remained  fixed  during 
each  run.  The  sinusoidal  acceleration  amplitude  was  increased 
at  a  slow  rate  from  a  static  position  to  vibration  tolerance. 

The  tolerance  criteria  was  aased  upon  individual  assessment 
concerning  a  combination  of;  vibration  severity,  pain,  loss  of 
physical  control,  or  blurred  vision.  When  tolerance  was  reached, 
the  subject  actuated  a  buzzer  which  began  the  RMS  data  collection 


procedure.  Each  subject  was  held  at  the  tolerance  level  for  20  to 
60  seconds.  The  time  interval  was  dependent  upon  the  stability  of 
RMS  data. 

To  gather  response  data  concerning  random  vibration, 
experiments  were  conducted  using  uniform  spectrum,  white  noise 
vibration  filtered  through  2  CPS  and  10  CPS,  bandwidth  filters. 

Random  vertical  and  angular  motion  data  were  recorded  and  are 
reported  in  Figures  2  through  4. 

The  tolerance  tests  established  the  upper  boundary 
of  human  response  in  accordance  with  the  experimental  control 
described.  To  validate  the  second  hypothesis  that  whole  body  response 
is  linear  below  this  level,  tests  were  performed  at  various 
acceleration  levels  down  to  very  low  intensity.  Regression  analysis 
indicated  strong  linearity  throughout  the  frequency  spectrum. 

Later  experiments  approached  'linearity  using  random 
waveforms  in  a  manner  analogous  to  that  of  a  constant  rate  spring. 
These  experiments  and  others  were  sufficiently  rigorous  that  within 
the  reasonable  constraints  of  our  test,  human  response  was 
sufficiently  linear  to  determine  the  transfer  function.  The  data  used 
were  mean  values  of  pooled  information  collected  at  the  seven 
different  acceleration  levels  (0.35  to  2.0  g‘s  P-P). 

Figure  5  describes  the  t ransmi ssibi 1  i ty  plot  for  mean 
values  and  the  90%  confidence  interval  of  head  acceleration  to  the 
input  acceleration. 


The  development  of  the  analytical  transfer  function 


was  achieved  using  Figure  5  and  an  asymptotic  approximation  method 

2 

of  the  frequency  spectrum. 

Equation  7  is  the  transfer  function  for  the 
acceleration  ratio  described  in  Equation  6  where  Cj  through  Cj2  are 
the  constant  coefficients  listed  in  Table  1. 
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CA(S)  =  TT  0j(.S)  =^£ut  {s) 

input 

i  =  1 

c,s5-  C2s4  4-  C3S3  ♦  C4S2  +  C£  +  C6 

s6  ♦  c7sJ  ♦  C8s4  ♦V-.+  c10»,*c1|»  +  c11 
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TABLE  I 


CONSTANT  VALUE 

Cj  n.28  X  JO1 

C2  10.40  x  104 

C3  38.70  x  '106 
C4  *  16.08  x  !08 

C5  44.31  x  1011 

C6  \  98.41  x  10H 


CONSTANT  VALUE 

C7  95.24  x  10l 
Cg  '  46.20  X  104 
C9  45.12  x  106 
C,0  10.17  x109 
C,,  64.20  x  1011 
C)2  98.42  x  10H 


Satisfactory  correlation  of  Equation  7  comparing 
analytical  results  with  experimental  measurements  was  achieved  for 
random  vibration,  sinusoidal  motion,  squared  shape  inputs,  and 
triangular  waves.  Correlation  was  secured  In  the  time  domain  and 
for  statistical  samples  (10  men}. 


•i’TTTrrrrr.m^n 


Frequently,  scaling  of  sensation  or  subjective 
response  is  approached  by  the  method  of  fractionization.  The 
stimulus  is  increased  or  decreased,  and  the  subject  reports  the 
point  where  the  sensation  appears  to  halve,  double,  etc.  In 
vibration  environments,  halving  or  doubling  the  stimulus  does  not 
halve  or  double  the  sensation.  In  human  vibration,  the  factors 
which  regulate  sensation  and  which  can  be  treated  analytically  have 
not  been  identified. 

It  was  our  third  hypothesis  that  under  vibration- the 
human  body  may  behave  in  an  elastic  or  inelastic  fashion.  When 
the  vibratory  energy  distorts  the  elastic  body,  dimensional  changes 
take  place  producing  reactions  which  tend  to  restore  the  body  to  the 
original  position.  The  work  that  is  performed  in  this  process 
balances  the  applied  load.  Consequently,  the  body’s  elasticity 
produces  restoring  forces  which  are  related  to  displacement.  The 
body's  vibratory  motion  continues  until  the  energy  imparted  is 
dissipated  or  removed.  The  time  rate  at  which  this  energy  is  used 
shall  be  referred  to  as  "absorbed  power".  "Absorbed  power"  is  a 
parameter  relating  vibratory  input  conditions  to  subjective  response 
and  is  defined  in  Equation  1. 

"Absorbed  power"  may  be  easily  produced  experimentally 
as  follows:  A  time  history  signal  of  acceleration  from  magnetic 
tape  or  directly  from  an  accelerometer  (during  field  measurements) 
is  sent  into  an  electric  circuit  representing  the  product  of  the 
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G.  (S)  and  the  Gp.(S)  transfer  functions.  The  output  from  this 
circuit  is  force.  The  input  acceleration  record  is  also  sent 
through  an  integrator  producing  velocity.  The  force  and. velocity 
signals  are  then  multiplied  electronically  and  put  through  an 
averaging  circuit  producing  average  "absorbed  power".  In  pre¬ 
liminary  experiments,  it  was  verified  that  "absorbed  power"  was 
zero  for  artificial,  inelastic  systems,  and  changed  exponentially 
for  characteri sti cs  representing  human  response. 

It  was  the  objective  in  the  initial  tolerance 
tests  (Figure  1)  to  secure  subjective  judgement  based  upon  a 
uniform  criteria.  From  these  tests,  mean  values  of  acceleration 
at  corresponding  frequencies  were  used  to  develop  a  constant  power 
curve  of  5.3  watts.  The  electrical  unit  of  "watts"  was  selected 
for  initial  use  in  that  it  was  readily  measured  and  described  in 
comparison  with  other  units  of  measure  (HP  and  BTU).  This  curve 
was  then  fitted  to  the  sinusoidal  tolerance  data.  Results 
indicated  that  "absorbed  power"  displayed  the  same  general 
characteristics  as  the  judgement  criteria. 

The  concept  of  "absorbed  power"  was  then  tested 
for  sensitivity  and  variability.  Experiments  were  designed  to 
evaluate  the  ability  to  discriminate  different  vibration  inputs 
and  to  order  their  comfort  or  severity. 

The  test  that  is  reported  here  was  performed 
using  twenty-six  random  waveforms  called  "rides".  Each  ride 


contained  different  PSD  vibration  characteristics  in  the  frequency 
spectrum  of  k  to  36  CPS  at  different  levels  of  acceleration.  The 
average  acceleration  used  to  order  the  rides  was  initially 
computed  using  four  control  subjects.  The  rides  were  presented  in 
sequence  of  increasing  acceleration  magnitude.  In  essence,  the 
first  ride  was  empirically  believed  to  be  the  smoothest  and  most 
comfortable  and  the  last  ride,  the  roughest.  A  group  of  ten 
subjects  judged  each  ride  (one  minute  duration)  subjectively  for 
comfort  and  reported  only  those  rides  that  appeared  to  be  "out  of 
place".  Acceleration  (g's)  and  power  (watts)  were  measured  during 
the  experiment. 

The  ten  subjects  reported  disagreement  in  60 
instances  with  the  acceleration  ride  sequence.  By  setting  three  or 
more  reports  per  ride  as  an  indication  of  significance,  nine  rides 
were  repeatedly  identified  as  misplaced.  Two  viewpoints  developed: 
(1)  either  the  rides  were  ordered  incorrectly,  or  (2)  the  subjects 
were  in  error.  The  "absorbed  power"  data  which  was  experimental ly 
measured  was  then  processed  and  compared  to  acceleration  (Table  2). 
"Absorbed  power"  closely  agreed  with  subjective  response  indicating 
the  acceleration  sequence  was  incorrect.  In  all  but  one  case  (Ride 
No.  10),  the  "absorbed  power"  measurement  appeared  to  identify  and 
discriminate  ride  comfort  for  small  as  well  as  large  differences. 

It  was  observed  that  the  measured  input  acceleration  did  not 
increase  in  precisely  the  same  ascending  order.  This  small 
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deviation  was  attributed  to  the  increase  in  the  test  population  from 
the  original  four  control  subjects  to  the  final  ten  subjects. 

\  TABLE  2 

Acceleration  -  "Absorbed  Power"  Comparison 


Ride  No. 

Input 

Measured 

No.  of 

Power 

( R.M.S.q's) 

Accel  oration 

Calls 

(Watts) 

Y 

0.417 

0.430 

\ 

0 

0 . 1964 

2 

0.463 

0.475 

\  0 

0 . 6864 

3 

0.479 

0.480 

\  .5 

.  V  •  ■  0 

0 .1742 

4 

0.506 

0.512 

0.3183 

5 

0.555 

0.559 

4 

0  .1635 

6 

0.6  58 

0.675 

0  • 

2.0976 

7 

0 . 660 

0.672 

\  6 

0 .3098 

8 

■  0.685 

0.703  ’ 

0 

0 . 5649 

9 

0.812 

0.828- 

7 

0.2657 

10 

0.840 

0 .850 

1 

0 .8080 

1 1 

0.841 

0 .881 

0 

3.6416 

12 

0.853 

0.855 

3 

2.2843 

13 

0.929 

0.924 

5 

0  . 6626 

14 

0.956 

'  0.962 

0 

1 .1632 

15  ■ 

1.007 

1.025  • 

2 

0 .6667 

16 

.  1.074 

1.121  ■ 

5 

0 .6203 

17 

1.144 

•  1.164 

0 

5.6228 

18 

1.162 

1.201 

5 

1 .6706 

19 

1.194 

1.205 

0 

1 .7719 

•  20 

1.205 

1  .21  1 

4 

0  .5415 

21 

1 .208 

1  .210 

0 

0  .9220 

22  • 

1 .438 

1.504 

0 

1.2092 

23 

1.475 

1 .628- 

10 

0  .8324 

24 

1 .617 

1.599 

0 

1.7269 

25 

1.808 

1.854 

0 

4.0612 

26 

1 .974 

1.942 

3 

1 .3736 

To  increase  the  utility  of  the  "absorbed  power" 
statement.  Equation  1,  research  work  was  directed  towards  developing 
a  simplified  relationship  which  would  not  requi re  stringent  control 
of.  the  input  data  and  which  also  would  include  the  effects  of  "time". 

It  was  rationalized  that  either  RMS  of  acceleration 
or  power  spectral  density  would  handle  all  the  input  requirements 
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which  may  be  defined  in  the  time  domain  or  statistically.  Equation 

8  is  the  new  expression  for  "absorbed  power".  Detailed  development 

is  given  in  Appendix  B.  It  describes  the  important  relationship 

that  the  total  "absorbed  power"  for  a  vibration  spectrum  is  the 

summation  of  the  power  at  each  frequency.  The  power  is  computed  as 

,  2 

the  product  of  the  acceleration  squared  (A.  r.m.s.)  at  each 
frequency  by  the  new  parameter  (K, ).  This  parameter  is  a  function 
of  frequency  but  does  not  vary  at  any  one  frequency. 


P(ave) 


„  .2 

K.A.  r.m.s. 
i  i 


Equation  8  describes  power  for  short  duration 
exposures  where  "absorbed  power"  is  not  affected  by  fatigue  or 
time  differences.  Time  alters  "absorbed  power"  and  this  effect 
is  included  in  Equation  9, 


(9) 


P  is  equal  to  the  sum  of  the  short  term  and  time 
dependent  "absorbed  power".  The  time  factor  l/t  represents 
the  approximate  onset  of  fatigue  and  t  is  the  time  of  the  exposure. 
Power  need  not  remain  constant.  It  may  be  a  function  of  time.  , 
However,  where  power  is  constant,  the  acceleration  relationship 
is  given  ass 


A 


(10) 
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Summary  and  Conclusions!  The  basic  premise  of  this 

i 

research  program  has  been  to  increase  the  knowledge  and  understanding 

of  human  response  to  mechanical  vibration*  Initial  research  steps 

were  basic  and  conventional,  and  centered  Upon  vertical  sinusoidal  ■ 

•  •  *  • 
tolerance  studies.  ! 


Experiments  followed  to  obtain  tolerance  data  for 
random  vibration  in  the  vertical  and  the  angular  modes  of  pitch  and 
roll.  After  this  accomplishment,  a  temporary  impasse  was  reached. 
The  problem  arose  of  describing  random  vibration  environments  which 
may  exist  in  an  infinite  number  of  patterns.,  PSD  tpower  Spectral 

i  \  •* 

density)  techniques  appeared  to  be  of  minimal  help.  A  new  approach 

was  selected.  The  validity  and  applicability  df  transfer  function 

'  .  '  ,  *  '  '  '  ' 

,  t  » 

techniques  to  human  dynamics  was  tested. 

Success  in  relating  analytical  transfer  function 

»  *  \  ■  ' ;  , ) 

,  .{  /** 

statements  and  experimental  test  data  was  reported  in  May  1965*  . 

Correlation  was  shown  for  white  noise  (random  vibration)  using 

W*  •*  i 

2  CPS  and  10  CPS  bandwidth  filters.  '  ’ 

i  The  value  of  this  analytical  approach  is  clearly 
demonstrated  in  synthesizing  such  items  as  effective  mass,  v  * 
impedance,  ^spring  rate,  etc.  from  two  basic  expressions  for  ,  ,t 

acceleration^  and  force*  A  case  in  point  is  the  impedance  transfer 

function  which  may  be  derived  as  fol  lows!  5  ,  ,  1  :  ,V 

.  ..  :Vi : .  V < 

•  v' '  v.;  •' •  V  ■ '  ‘  '•'  'l 


■it'  •;  .i  .•••  • 


v:  : 

I  '•  1  '  ’  <  *  ’ 


•i 

X.  ’:  ::v  {■; 

» I  *  •i' 


.•  -  ;  /  I'.'.:'::-; 

^  .'-iT  ;  V.;;vf  ■  •; 


V'm  • 


i  ■  v  .  /»  *  • 

. '  <■  i  -i  •  . 

'  ‘ ■  >  • .  ■  • 

■  i  r 

}■  ■  i  •  )  «• 

;  •••  •  •; 

•  j-  ■„  .'•*  ,;v 

•  Si  ••  •. 
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IMPEDANCE 

GZ(S) 

= 

input 

Si  nee 

A(S)  . 

=  SV(S), 

and  acceleration 

GA(S) 

=  *head-(s) 

input 

Force 

Gp(S) 

=  ^3=(s) 

head 

GZ(S) 

=  SGf(S)  Ga(S) 

Expanding  the  transfer  function  approach,  the  new 
criteria  "absorbed  power"  was  developed  by  Mr.  R.  A.  Lee.  It  was 
shown  that  "absorbed  power"  corresponded  to  subjective  response. 

Later,  it  was  clearly  indicated  that  subj  ecti  ve '.response  can  be 

\ 

forecasted  using  "absorbed  power"  and  that  the  effect  of  time  may 
be  treated  analytically.  \ 

Detailed  examination  of  the  "absorbed  power" 
criteria  indicated  advantages  not  present  in  the  "acceleration" 
measurement.  Acceleration  is  a  vector  quantity.  It  must  be  assessed 
in  combination  with  frequency  spectrum,  magnitude,  and  direction.  . 
"Absorbed  power",  however,  is  a  scalar  quantity  which  may  be 
described  by  magnitude  only.  Thus  it  is  additive  and  may  be 
summed -in  mul ti -deg ree -of -freedom  environments. 


Acceleration  data  measurements  are  usually  collected 
at  input  locations  where  the  influence  of  individual  differences, 
posture,  or  seat  cushions  cannot  be  separately  assessed.  In  this 
important  area  where  acceleration  demonstrates  insensitivity, 
"absorbed  power"  displays  sensitivity  and  varies  with  experimental 
and  anatomical  differences. 

Figure  6  shows  "absorbed  power"  versus  ride.  The 
thirteen  rides  plotted  on  the  abscissa  are  random  records.  Each 
ride  represents  different  descriptions  of  energy  within  the  same 
frequency  spectrum.  It  can  be  seen  that  the  power  which  was 
analytically  determined  corresponds  closely  to  the  experimentally 
measured  "absorbed  power".  Figure  7  describes  the  application 
of  Equation  10  and  compares  the  reciprocal  of  "absorbed  power" 
(increased  magnitude  decreased  severity)  to  the  time  of  exposure. 
Figure  7  reveals  the  utility  of  this  analytical  approach  to  define 
time  of  exposure  as  well  as  ride  severity. 

Thus,  the  following  distinct  and  separate  conclusions 

are  drawn: 

1.  Human  response  to  sinusoidal  and  random  vibration  displays 
linear  characteristics  under  physical  equilibrium. 

2.  Transfer  function  statements  accurately  describe  human  response 
to  random  vibration. 

3.  Human  vibration  may  be  studied  analytically  using  transfer 


function  statements. 
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“Absorbed  power"  describes  human  response  quantitatively  and 
is  sensitive  to  time. 

Experimental  determination  of  “absorbed  power11- does  not  require 
frequency  spectrum  analysis. 

Optimization  studies  involving  human  dynamics  may  be  conducted 

\ 

continuously  without, time  lapse  that  occurs  for  frequency 
spectrum  analysis,  \ 
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Appendix  A 


TABLE  I  -  Statistical  Data 


Standard 


Averaqe 

Devi  ati  on 

Age  (years) 

'  28.00 

10.36 

Height  (inches) 

Height,  stature 

70.27 

2.64 

1A-  Sitting  height,  erect 

35.70 

\ 

1.92 

IB.  Eye  height,  normal  sitting, 

30.90 

2.05 

canthus 

1C.  Buttock-shoulder  height 

23.81 

1 .76 

ID.  Buttock-elbow  height 

9.76 

1.62 

IE.  Seat -height 

16.69 

0.98 

IF.  Knee  height 

22.66 

0.98 

Weight  -  no  equipment  (pounds) 

175.17 

22.19 

1.  Sitting  weight,  erect 

148.78 

11.19 

2.  Sitting  weight,  braced 

114.71 

21.54 

Trunk  (inches) 

2A.  Shoulder  width 

18.33 

0.93 

28.  Elbow  wi  dth 

19.47 

1.50 

2C.  Seat  width 

14.86 

1.28 

3A-  Chest  depth 

9.76 

1.18 

38.  Abdomina.1  depth 

00 

• 

cr\ 

0.81 

Head  (inches) 

Head  length,  front  to  back 

7.73 

0.29 

Head  width,  side  to  side 

6.26 

0.21 

Hand  (inches) 

7A  Hand  length 

7.93 

0.51 

O'? 


7B 

Hand  width 

3.69 

0.31 

Arm  (in 

ches) 

4A 

Elbow  -  finger  length 

19.34 

0.92 

A -2 
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APPENDIX  B  Derivation  of  Equations 


The  equation  for  RMS  of  acceleration  is: 


rms 


1  im  _J_ 

T 

!  T  -* 


QO 


/ 


'(t)  dt 


where  a(t)  is  any  acceleration  as  a  function  of  time  and  A  is 
1  rms 

the  RMS  value  of  a(t). 

If  displacement  is  sinusoidal  and  is  given  by 
_  Ap 

s - ^  sin  cot 

co 


Then  velocity  is  given  by 


^  =  V(t)  =  cos  co  t 


A 


where  A  is  the  peak  amplitude  of  the  accelerations  (from  zero  to 
P 

peak)  then  the  acceleration  for  this  velocity  is 
~  =  a(t)  =  -A  sin  cot. 


and  the  RMS  of  this  sinusoidal  acceleration  from  (l)  and  (3)  is 


A  =  Hm  1 

rms  ->■  oo  — 

Aj  T 

Squaring  both  sides 

,2 


A^  sin^cot  dt 


A 


rms 


2 
A 

OO  -y 


1  i  m  A 

T  -*■  —  -E. 


I  sin' 


'cot  dt 


Center 


0) 


(2) 


(3) 


(4) 


(3 


o 


B-3 


Then : 


P  = 


lim  J_ 

^■oo  T 


Now 


Va2 


(sinco't  CO  Slot  cos  f  +  cos2 


<ot  si 


n  if) 


dt 


J~  Si"  “t  dt  =  -  i”-"  aot  and/cos2Ut  dt  =  it  +  ^  sin2ut 


(10) 


Therefo re 
1  im 


T 


and  since 


1 1  m 


T 


03 


1  T 


cq 


■I 


sin2ut  dt  =  _L 

T  -*■  oo  T 


T 


cos  cot  dt 


i 

f  sirwt  cos-ut  =  _llm  i 

K  T  T 


f  II  j  S1’n2wt  dt  = 


lim  ]  r  ]  T 

oo  T  L"I^  cos2cot]  =  0 


01) 


(12) 


the  equation  for  average  power  becomes 
P  =  G  si  n  if 


co 


lim  A 

T  ->  — £  r 

co  j  \J 

P 


2 

si  n  cot 


Equation  M  E^tf°"  “3)  «» 


same  as 


rms 


Therefore,  Equation  (13)  can  be  wri 


wn  tten 


P  = 


Kp  si 


i  n  if 


A 


(0 


rms 


(13) 


(14) 


Since  Kq  and  f)  are  constant  at  one  frequency.  Equation 
•A;  can  be  written 


B-3 


Using  the  identities 

i  r 

sin  at  cos  bt  =  j  sin  (a-b)t  +  sin  (a+b)t 


cos  at  cos  bt  =  y 


cos  (a-b)t  +  cos  (a+b)t  J 


one  obtai ns 


1  i  m  1 


oo 


f  /  si 


sin  at  cos  bt  dt 


1  i  m  j_  cos(  a-b)t  cos(  a+b)t 

T  oo  T  L  2  (a-b)  2(a  +  b) 


=  0 


1  i  m 

T  ->  oo  T 


1  i  m  J_ 
T-o)  T 


cos  at  COS  bt  dt 


sin( a-b)t  +  sin  (a+b)t 
2  (a-b)  2( a+b) 


=  0  a  /  b 


(24) 


(25) 


and  since 
1  i  m 


"1"  oo  T 


cos  oo.  t  cos  co.  t  dt 
1  J 


lim  i  r  2  . 

T  —1 — /  cos  oj.t  dt:  i  =  i 

T  oo  T  j 


(26) 


and  Equation  (12)  holds.  Equation  (23)  becomes 
n 


'  P  = 


F.V.  sin 

ii  /  i 


i  = 


r  i  ,•  ~  i  t  si  n  2  co.T 
I  nm  1  /T  i 

T  2  TST - 


) 


i  = 


^i ^i  sin  L P  . 


ir  1  i  m 

1  1 

/  .  ^  \  I 

oo 

2  +  T 

y  si  n  2co^T  J 

4co. 


(27) 


b-6 


'-/V 


0 


Taking  the  limits  one  obtains 


T  V  F.V.  sin  Ip . 
2  ii  Vi 


i  - 


Since' for  sinusoidal  waves 


F.  =  A  V.  =  A  /co.  and  A2  =  2A2  * 

”  ~  “  p  rms 


i  Gi  Pi  i  P/  i 


and  using  equation  (16), 

Equation(28)  can  be  written 


■r 


k.  k: 

i  i 


rms 


i  =  i 


F  rom  ( 5) 

,2 


i/ 


*  A“  1  i  m 

rms  =  i  _>  f  J  P 

05  T  o 

rA2  t  A 

_E 


2  2 
A  sin  cot  dt 


A 


’  _  1  i  m  ]_ 

rms  T  ->50  T  L  2 


2  T 

kf-  sin  2ut 


A 


1  i  m 


rms  T 


GO 


r  2  2 

I  A  -  A  si  n  2  coj 

L  — E  _E 

2  Zrot 


A*  -  1e  or  A2  =  2A2 
rms  — c-  r> 


rms 


(28) 


(29) 


(30) 
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If  the  exponent  is  odd,  the  form  is 


R  *  l-D  '2'  IkV*-1 

.  .  ...  i_ 


*  r~3  .  v  r-5  -  +  .  .  .  (35) 

K,u  +  r  f 

O  P  •  8 


=  (-] )  2.  W  ju 


+  K^co1"""^  -  +  .  .  . 


Where  r  is  the  degree  of  the  polynomial  in  either  the 
numerator  or  denominator. 

Writing  Equation  (32)  in  the  form  - 

K  (F.  +  JcoFp) 

«•“.»  -  \  I  s 


where  the  F's  are  given  by  Equations  (33)  thru  (3&) 

Multiplying  denominator  and  numerator  of  Equation  (37)  by 

(f^  -  JuF^)  yields  _ j 

a(Jw)  =  K  !  F1F3  *  ^  F2F4  *  (F2F3  "  F1F4}  ;  (38) 

0  i  o  2  2  '  | 

|  F3/  h-  (/  FkZ  _J 


where  x  =  K 


G(Ja>)  may  be  thought  of  as  being  in  the  form 


G(Ju)  =  x  +  yi 
2 

F.F0  +  to  F^Fi 

,  t  3 _ £_Ji 

0  .  2  2^  2 
+  “  f4 


y  =  Kc 


<'V3  -  fifa) 

F32  +  u2F42 


The  number  V=  x  +  y  is  called  the  modulus  or  absolute  value 
of  the  complex  number.  The  smallest  positive  angle  which  the  vector  makes 
with  the  positive  x-axis  i  s  -.represented  by  Arctan 

x 


Hence 


G(  Ju) 


!Ko2  (F,F,  ^V4)2.,,  2  PH  F2F,  -  F,F4)2 

A  1  (c  2  ,  2  2,2  °,c  2  2  2,2 

\j  (F^  +  w  )  ^3  +  u  F4  ) 


lG  (Ju) 


Kq  .  |  (F]F3  +  u2F2F4)2  +  u2(F2F3-F]F4)2 


7\  j 

i  ■ 

V 


,,  2  2  2,2 

(r3  +  u  F4  ) 


j G( Ju) 


=  K, 


!F3  +  ^2F/4)2  +  u2  ( F2F3  -  f1f4)2 
(F32  +  u2F42) 


(39) 


The  angle  which  the  vector  makes  with  the  x-axis  is 


vj?  .  ,  1  ^F2F3  "  F1*V 

;  —  Arctan  . - x — 


/  ' 


L 


F1F3  +w  F2F4i 


Equation  (39)  may  be  simplified  as  follows 


G(  Ju) 


.  1  2  2  2  4222222  222 

i  ,  \  F1  F3  +2“  F1F2F3F4^  F2  F4  ^  F2  A  -2“  FlF2F3F4+a)  F]  F4 

1  7  F32  +  “V 
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J  r  2  2  2r  2 *  2  2,  2_  2  _  2. 

\  F ]  (F3  +  )  +  u  F2  (fa)  l~4  +  F3  ) 


|uWcu;i  2  2  2 

F3  +  “  f4 

I  _  2““  ™  2“~2-' 

j  ! _  YF3  +  w  t-,.  )  ( F 1  +  u  r2  ) 


>  G  (  Jco ) 
i 


2  2  2 
F0  +  u  F, 

3  4 


:  2  2_  2 
F1  +("F2 
?  c  2  .  2  2 
j  F3  +  03  f4 


If  G(Ju)  is  a  transfer  function  defined  by 

Force  _  , 

Acceleration  '  u 


F.  =  A.  G(  Ju.  ) 

i  i  i 

rms  rms 

Using  Equation  (28)  and  the  relationship  given  in 
Equation  (29),  ‘'absorbed  power"  can  be  written  as 


l  l  ^ 

P  =  V  1  ( F.A  /  )  sin  J  i 

4 ,  2  1  1 


Substituting  Equation  (42)  into  Equation  (43), 


n  |  |g(  J^)[  A  2 
P  -V*  _ rms 

L  L  «• 

i  =  1 


.  ip  , 

sin  i 


This  is  the  same  as  Equation  (30)  with 


K.  = 

i 


•  G  ( Jco .  .j 
i  1 


B  —  1 0 


r 


«  <  i  * 


where 


G(  Jw)  j  ~  K< 


A  S 


!  _  2  2  2 
i  +  a)  F2 


2  2  2 
\|  F3  +  “  F4 


(46) 


And 


.  (/)  w(F1F4“F2F3') 

i  n  =  - - — 


si  n 


3  •  “  ‘  2‘  4 

Squaring  and  cancelling  yields 


Aj  (F j F,  +co2F0F;,)2  +  uAf^.-F^F,)2 


1‘4  ‘2*  3' 


st  n 


=  m(F1F4~F2F3^ 


L  2,  2  2r  2.  22,  2_  2  _  2. 

't*l  (F^  +  0  F^  )  +  «  F2  (10  +  F^  ) 


si  n  GA 


^F,VF2F3i 


/\  ,  2_22,  2  ,  2C  2. 

/\!(w  F^  +  F3  )  ^ F 1  +  u  F2  ) 

Substituting  the  relationship  for  jG(Jco)  |  and  sin 
yields  ^ 


K.  = 


F,  F,  -  F„F. 
i  4  /  j 

f32  +  "V  ! 


(47) 

into  Equation  (45) 

(48) 


Using  Equation  (48)  and  equation  (30),  an  "absorbed  power" 
can  be  calculated  for  any  given  frequency  spectrum. 


-I 


F ,  F, 


F2F3 


A. 


(49) 


i  =1 


2_  2 
+  0)  1- , 

H 


rms 


B-l  1 


J 


Solving  for  A.  in  Equation  (30)  gives 


or  any  given  power  and  frequency,  an  acceleration  can  be 
ca i cul a ted. 


(50) 


For  vibrations  of  short  duration,  “absorbed  power"  can  be 
used  as  a  measurement  of  vibration  severity.  To  include  fatigue, 
the  energy  absorbed  is  modified  to  include  the  time  at  which 
fatigue  begins. 


=  P  + 


Uo 


i  ?  dt 


(51) 


ere  is  the  long-term  "absorbed  power",  p  is  the  average  pov/r:-. 


—  is  the  time  scale  factor 

t 


the  initial  time  t  is  the  onset  of 

o 


fatigue.  Time  t  is  the  length  of  exposure. 

Power  need  not  remain  constant  but  can  be  a  function  of 
time.  If  power  remains  constant,  wa  obtain 


?  a  ?  p  - 

L  L 


(52) 


Substituting  KA  for  power. 


2  2  t 

?.  =  KA  +  KA  —4 

t  t 


(53) 


And  an  acceleration-time  relationship  is  given  by 


o 


(54) 


